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ABSTRACT 

It is pointed out that electrical surge-generator circuits are, essentially, periodi- 
cally loaded lines. It is then shown that, with the aid of certain reasonable 
assumptions, the well-known theory dealing with such lines leads directly to 
expressions for the free periods and damping factors of surge-generator circuits. 
The resulting equations can be solved exactly in some cases and can be solved 
approximately in all cases. In certain cases exact expressions can be obtained 
for the current and voltage amplitudes. 
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I. INTRODUCTION 

Kecent investigations into the behavior of transmission lines and 
associated equipment, when subjected to voltage and current surges 
simulating those produced by lightning disturbances, have led to the 
development and use of the so-called surge-voltage and surge-current 
generators. A surge-voltage generator consists of a number of capa- 
citors which are charged in parallel and then, by a suitable arrange- 
ment of spark-gaps, discharged in series. Figure 1 is a schematic 
diagram of the circuits of a surge-generator as they function on dis- 
charge, the charging circuits being omitted to avoid confusing the 
figure. If all the capacitors are identical and are initially charged to 
the same voltage E , the voltage appearing across the terminals of 
the nth mesh will rise to a maximum value of the order of nE when 
the gaps are discharged in succession, beginning at the left-hand end. 
A surge-current generator consists of a number of capacitors which 
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are arranged to be charged and discharged in parallel. The current 
in the terminal mesh of an ^-section generator will be roughly n times 
that obtainable from a single capacitor. Figure 3 shows a schematic 
diagram of a surge-current generator in which, again, the charging 
circuit is not shown. 

The methods of obtaining approximate solutions for such circuits 
are obvious. If one neglects the internal structure of such a circuit 
and considers it to consist simply of a lumped effective capacitance, 
inductance, and resistance, the frequency and amplitude of the 
fundamental component of current or voltage can be computed with 
a reasonable degree of accuracy. As an additional approximation 
one might consider the entire generator to consist of two circuits, a 
generator with lumped constants coupled to a load circuit. Appar- 
ently no analyses of the discharge characteristics of surge-generator 
circuits have been published which have gone be3 r ond these two 
obvious approximations, although R. Eisner and C. S. Eoys have 
published complete analyses of the charging characteristics of such 
circuits. 1 While serving fairly well to predict the frequency of the 
fundamental component of current or voltage, this approximation 
leads to erroneous results when an attempt is made to use it in the 
computation of critical damping resistances, as will be shown in the 
section on surge-current generators. In addition, such a simplifica- 
tion entirely neglects the presence of frequencies higher than the 
fundamental. It is just these higher frequencies, however, which are 
of particular importance in defining the wave front of such surges. 

It is the purpose of this paper to point out that surge-generator 
circuits are essentially periodically loaded lines, and that the highly 
developed theory of the periodically loaded line can be applied directly 
to surge-generator problems. It is shown that straightforward calcu- 
lations lead to expressions for the free periods and damping factors of 
surge-generator circuits which can be solved exactly in some cases and 
approximately in all cases. In certain cases exact expressions can be 
obtained for the current or voltage amplitudes. 

II. SURGE- VOLTAGE GENERATOR 

1. STATEMENT OF THE GENERAL PROBLEM 

The problem of the surge-voltage generator can be stated conven- 
iently in terms of the simplified circuit shown in figure 1. A set of 
n-1 identical capacitors of capacitance C are first charged in parallel to 
a common voltage E . By a suitable arrangement of sphere gaps, 
as shown, they are then discharged in series through leads having 
identical resistances R and inductances L. The nth. mesh, of resist- 
ance R n , inductance I n , and capacitance C n , contains the equipment 
under test as well as other measuring apparatus, and is in general 
different from the other n-1 meshes. The capacitance to ground 
associated with each capacitor and its two spherical terminals is 
supposed, for simplicity, to be concentrated as shown in figure 1 and 
to be of common value G g . Different locations of this lumped 

i P. L. Bellaschi, Trans.'Am. Inst. Elec. Engrs. 51, 936-945 (1932); 53, 544-552 (1933); Elec. Eng. 53, 
86-94 (1934). $ to 

J. C. Dowell and C. M. Foust, Trans. Am. Inst. Elec. Engrs. 52, 537-543 (1933). 

B. L. Goodlet, J. Inst. Elec. Eng. (London), 74, 377-396 (1934). 

K. B. McEachron and J. L. Thomason, Gen. Elec. Rev. 38, 126-131 (1935). 
R. Eisner, Arch. Elektrotech. 29, 655-682 (1935). 

C. S. Roys, Engineering Bulletin Purdue University 19, no. 3 (May 1935). Research Series no. 50. 
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capacitance to ground, as for example, splitting it between the two 
capacitor terminals, would present no essential mathematical diffi- 
culties. We neglect the internal structure of the capacitors, if any,, 
and assume that their internal losses are sufficiently represented by 
a suitable increase in the lead resistance R. 

All the capacitors are initially charged to the same voltage, E Q . If 
now the first gap is broken down by some artificial means, each of 
the other gaps in the circuit will break down in its turn as the voltage 
across it rises to its breakdown value. The problem, then, is to find 
the charges on the various capacitors as explicit functions of time, 
subject to the given initial conditions of charge and current. Let 
Q s represent the charge on the sth main capacitor and q s the charge 
on the sth capacitor to ground. The method of solution is to set up 
the equation for the summation of the voltages around each mesh in 
terms of the charges on the various capacitors. The q's can be elimi- 
nated in terms of the Q's. After certain transformations the result- 
ing n second-order differential equations in the Q's lead to a trigono- 
metric equation of degree n, the secular equation, the roots of which 
^furnish the free periods and damping factors of the system. These 
roots can be obtained exactly in some cases and can be approximated 
in every case. The 2n constants necessary to satisfy the original n 
second-order differential equations determine the amplitudes and 
phase angles of the charge components. Exact expressions can be 
found for these constants in certain limiting cases and approximate 
expressions in certain others. 

2. FORMAL SOLUTION OF THE GENERAL PROBLEM 

Keferring now to figure 1 and writing down the equations for the 
summation of the voltages around each mesh, the gaps now being 
considered conducting, 
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Figure 1. — Schematic diagram of a surge-voltage generator. 
The charging circuits have been omitted to avoid confusing the figure. 



-Bi l -Ldii/dt+Q 1 /O+q x IC t =0'\ 
-blCi-Rh-Idk/dt+Qj/C+tolC^O 



-q,- 1 /C-Ri,-Ldi,fdt+Q,/C+q./C,=o\ 
— g»_i/C— R»% — L n di„Jdt + QJC n = 0. 



(1) 



Now consider the currents at a representative point, as a. By definition 

h=-dQJdt (2) 
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Then at point a 

—dQ s /dt+d(i s /di= —dQ s+i fdt 



Integrating, 
At *=0, 



q.s=Q s —Qs+i+h s (constant) 



Qs—Qs,o 
Qs+i—Qs+1,0 



(3) 



Therefore, 

At=2*,o — Qs,o~rQs+i,o 
Substituting eq 2 and 3 in eq 1 we have, 

LdtQi/d?+RdQ 1 ldt+Q l IC+Q 1 /C,-Q l /C,+h 1 /C § =& 
- QilC,+Q a /C-h l lC t +Ld i Q a fdt i +BdQJdt+Q i fO 

+Q i {C,-Q t /C,+h i /C,=0 



(4) 



-Q.-ilC t +Q,/C,-h,- 1 IC,+L(PQ,/dP+RdQ,ldt+Q,/C 

+Q.fC,-Q, +1 /C,+h,/O,=0 



(5) 



-Qn-ilO e +Q n /C-h n ^IC +LJ 2 Q n /df+R n dQ n /dt+Q n /C a =0) 

Now make the usual exponential substitution 

Q,=d,+A*- M (6) 

where d s , A s , and X are constants, as yet undetermined. 2 Applying 
eq 6 to the typical equation of eq 5 we obtain, 

-h s - 1 -d s _ l +[-A s .. 1 +(C e L-K 2 -C,R\+CJC+2)A 1 - 
+ (C,/C+2)d,-d, +l +h.=0 

This suggests that we define two new quantities 



-A s+l ]e 



w(\)=C e L\ 2 -G e R\+C s /0+2 \ 
X(\)=C s L tt \ 2 -C ! fi n \+C l ,/C n + 1/ 

Eq 5 then become 

[(w-l)A 1 -A 2 ]e-^'=-(G e IG+l)d 1 +d 2 -h 1 
[-A l +wA 2 -A 3 ]e-»=d 1 -(CJC+2)d 2 +d z +h l -h i 



(7) 



[—A s - 1 +wA s —A s+l ]e- > ' t =d s - 1 — (C IC+2)d s +d 3+1 +h s - 1 — h s 



l-A.- 1 +XA n \e-"=d^ 1 -(C,IC+l)d n +h n - 1 



(8) 



Since the constants d are at our disposal, we shall so choose them 
as to reduce the right-hand members of eq 8 to zero. The efs are 
therefore defined by the following equations, 



271 

2 Strictly speaking there will be 2n \'s and Q, will be of the form Q s = ^2 A ^ e ~^ kt ' For the time being 

k = l 
we shall omit the subscript k. Strickly speaking also A, is not necessarily a constant. In case a given X r is 
a double root of the secular equation, that is, if the rth frequency component is critically damped, the 
corresponding A ir is a linear function of time of the form {At-\-B) . See, for example, Ames and Murnaghan, 
Theoretical Mechanics, page 124 (Ginn and Co., 1929). 
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(9) 



With the d's having the values determined by eq 9 the solution 
of eq 8 reduces to the solution of the following equations: 



(w-l)A 1 -A 2 =0) 

—A l +wA 2 —A z =0 



(10) 



—As-i+wAs— A s+1 =0 

' '-A n ^+XA n =0) 

The condition that a solution other than zero shall exist for the A's 
requires, as usual, that the determinant of the coefficients shall 
vanish. Let us call this determinant D n . The condition for a solution 
then reduces to 
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(11) 



(12) 



Expanding D n on its first column, we obtain 

D n =(w-l)H n ^-H n _ 2 (13) 

Expanding H n similarly, we obtain 

Hn^wH^-Hn^ (14) 

The method of handling such difference equations and the applica- 
tion of the results to electrical problems has been discussed by Wheeler 
and Murnaghan. 3 Let us make the trial substitution H n —cx n ) where 
c is any constant. Substituting in eq 14 we obtain 

x 2 — wx= — 1 
whence x = w/2 ± [w 2 /4 — 1 ] * /2 

This suggests the substitution 

w/2 = cosh (15) 

3 H. A. Wheeler and F. D. Murnaghan, Phil. Mag. [7] 6, 146-174 (July 1928). 
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Making this substitution, we obtain 

x = cosh 0±sinh 0=6=*=* 

Hence 

H n =cx n =ce :hn<i> 
or H n =CiCOsh ti0+c 2 sinh n<j> 

The values of c x and c 2 may be determined from the values assumed by 
H n when n=0 and n=l. From eq 12 and 14 we have H 2 =wX— 1 = 

Therefore 

flo=l 

Using these values we readily find 
X— cosh 



c 2 =- 



sinh 



so that 



And finally 



j-r X sinh n<j>— sinh (n— 1)0 /1c v 

ft " Sh0 (16) 



n _ Xcosh (n— 1/2)0— cosh (n— 3/2)0 ( . 

^ w ~ cosh 0/2 (17; 

The problem of solving eq 10 is therefore reduced to the problem 
of finding those values of which will reduce eq 17 to zero. There 
will be n such values of corresponding to the n degrees of freedom of 
the system. Exact solutions for will be found only for certain special 
values of the constants R ny L n , and O n) but approximate solutions 
can always be made. 

Values of X are obtained by substituting these values of in eq 7 
and (15), from which we obtain 

X w ,X* i2 =i2/2Z±[iP/4Z 2 ^ (18) 

there being n values of k corresponding to the n possible values of 0. 
It can now be verified by direct substitution that eq 10 will be 
satisfied for all possible values of if the A's are of the form 

-4,=P t cosh (s- l/2)fo (19) 

where the P k are constants. 

Therefore, the final solution of eq 5 is of the form 

Qs=d 8 +^[P ktl e-^ t +P kt2 e~^ t ] cosh (s- 1/2)0 (20) 

% 

The constant d s is evaluated from eq 9 and the constants P M and P k>2 
are determined by the initial conditions, that is, by the values of 
Q s and i s at 2=0. With the evaluation of these constants the problem 
is completely solved, for with the aid of eq 20 we can now compute the 
voltage drop across any element of the circuit of figure 1 as an explicit 
function of time. 
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3. ILLUSTRATIVE CASE, EXACT SOLUTION 

The application of the preceding formula is best illustrated by 
means of a specific example. For simplicity, a case has been chosen 
in which an exact solution is possible. Other special cases leading 
to exact solutions can be discovered with little difficulty. These 
exact solutions may be found to be of little practical importance since 
the constants of the terminal mesh, including as they do the constants 
of the apparatus under test, cannot in general be chosen arbitrarily. 

Referring to figure 1, let us assume that 

CG 
R n =K } ±j n =L) ^n— /i\p 

This is assuming that the terminal section is identical with the 
typical section, the capacitance C n being replaced by a capacitance 
C in series with a capacitance C g . Substituting these values in eq 7 
we find 

X=w=2 cosh 

Substituting this value in eq 17, we find 

n ^ cosh (n+l/2)0 
Un cosh cj>/2 

The condition for a solution is that D n =0, or 

cosh (n-hl/2)</>_ 



cosh <£/2 



= 



Hence ^ >k== 2n+l Kl 

or,if*0=0 0,=2^7r,/c=o,l,2,3, . . . n-1 (21) 

Referring to eq 18, we have 

KuK2==R/2L±[R 2 /4L 2 -l/CL-(±/C g L) sin 2 k /2]^ 2 (22) 

Three cases are to be distinguished, as usual, depending upon 
whether the radical of eq 22 is real, zero, or imaginary. In the first 
case (R greater than necessary for critical damping) the amplitude 
substitutions will be as given in eq 20. In the second case (R equal 
to the critical damping resistance for a given frequency) the sub- 
stitution will be of the form. (See footnote 2, page 588.) 

(P r ,lt + P rt2 )e~ X rt CQS (£-1/2)0,. 

In the third case (R less than that necessary for critical damping) the 
X Vwill be of the form 

X/ b =/i=b27r / i^ 

where \x—R\2L 

2wv!=[-R 2 /4L 2 +llCL+ (4/tf,L) sin %/2] 1 '* (23) 

In this case the substitution corresponding to eq 20 will be of the form 

Pfccos (s— lfflfar** cos (2Tv k t—a te ), 

where P k and a k are now the unknown constants. 
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Suppose now that the value of R is such as to critically damp the 
rth frequency component. We can then write a general expression 
for the charge on the sth capacitor in the following form: 

Q s =d s +£(P kfl e-^+P kt2 e-^) cos (s- 1/2)0* 

fc = 



n-1 

+ S Pjc COS (s— 1/2)<V~ M ' cos (2icv k t—a k ). 



(24) 



The 2n constants are determined from the initial conditions, that is 
at /=0, 

Qs=Q s ,o 

The problem of determining these 2n constants is the standard 

problem of passing an n term trigonometric sum through n arbitrarily 

chosen points. If the n points are equally spaced (that is, if the solu- 

kir 2&4-1 
tion of eq 17 yields exact values of 6 of the form 6= — > , . ir, etc.) 

n zn~\~ l 

the solutions can be written down immediately. 4 The condition that 
Q S =Q S at t=0 leads to the solutions of the following identical form 

/ 2&4-1 \ 

for I P*,i+P k , 2 ), P r ,2, and P k cos a k (for the case where fl t = 2 T. tJ 

P k cos a k =4/(2n+l)h(Qs,o-d s ) cos (*- 1/2)0*. (25) 

The second initial condition, i s =i s ,o at tf=0, leads to similar solutions 
for the same pairs of constants, from which the individual constants 
can be evaluated. If the initial conditions are known the summations 
of eq 25 can always be carried through and the constants determined. 
The exact determination of these initial conditions is a matter of 
some difficulty in the case of a surge-voltage generator. It is not 
correct to assume that at time zero all charges are equal and all 
currents zero, although this represents the limiting case of all gaps 
closing simultaneously. Actually, the breakdown of the gaps is 
successive, and the initial time for the whole system should be taken 
as the moment of breakdown of the last gap. At this moment the 
charges on the preceding capacitors are not identical and the currents 
in the preceding meshes are not zero. To compute the initial con- 
ditions exactly the following procedure must be followed. First, 
compute the voltage above ground of the terminal of the second 
section of the generator (that is, point a 7 fig. 1) after the first gap has 
been closed and before the second gap has broken down. This com- 
putation can be carried out exactly since there is only one gap in this 
circuit and it is arbitrarily closed at time zero, at which time both 
initial charges are equal and both initial currents are zero. Now, 
having a knowledge of the voltage at which the second gap will break 
down, or having made a reasonable assumption as to that voltage, we 
compute the charges and currents existing at the moment that break- 

* W. E. Byerly, Fourier's Series and Spherical Harmonics, chapter II. (Ginn and Co., 1895.) 
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down voltage is reached at the second gap. This moment is taken 
as the new origin of time measurement in dealing with the three- 
section network now in the circuit, and these charges and currents 
are the required initial conditions. We repeat the above process for 
the third gap and for each of the other gaps in turn until finally we 
are able to compute all the charges and all the currents at the moment 
of breakdown of the last gap. It is readily seen that the rigorous 
computation of the voltage across the terminals of an 7i-section gen- 
erator becomes very laborious if n is very large, say 20 or 30. 

To make these ideas somewhat more specific, sample calculations 
have been carried out on a three-section network having the following 
circuit constants: 



2.0 




MICRO SECONDS 

Figure 2. — Voltage surge of a two-section generator discharging into a previously 
uncharged terminal section. 

Curve I, both gaps closed simultaneously. 

Curve If, second gap closed 1 X 10 -7 seconds after first gap. 

L=L n =5 X 10- 6 h. 
R=E n =4:0.5 ohms. 

(7=0.25 X 10- 6 farads. 
C g =0.25 X 10- 8 farads. 

fl C Cg 

° n ~c+c g 

E = initial charging voltage. 

The value of resistance used is just sufficient to critically damp the 
fundamental of the three-section network. We have here a two- 
section generator discharging into a previously uncharged terminal 
circuit. 



594 Journal of Research oj the National Bureau of Standards [vol. n 

In figure 2 are plotted curves of E 3 /E against time, E 3 being the 
voltage across the capacitor C g in the terminal circuit. Curve i" is 
plotted for the limiting case in which both gaps are closed simultane- 
ously. Curve II is plotted for the case in which the second gap closes 
0.1 microsecond after the first gap. In this time, the voltage across 
the second gap will have risen to about 1.25 E . The difference be- 
tween these curves indicates the importance of initial conditions, that 
is, gap settings, on the shape of the voltage wave front. 

4. DISCUSSION 

In surge-voltage generators as used in practical work the situation 
is somewhat more complicated than is indicated in the simple example 
just cited. The inductance and capacitance of each section of an 
actual generator might be expected to be of the order of magnitude 
of those given in the example, but the resistance of each section, R, 
might be much less than that required for critical damping of the 
fundamental. An aperiodic wave is then obtained by inserting re- 
sistance in the terminal section, that is, making R n large. Ordinarily 
a resistance equal to several times the surge impedance of the circuit 
(Ltot&i/Ci) 1/2 , is found sufficient to give a reasonably smooth wave. 5 

A solution for the frequencies and damping factors is obtained, as 
before, through the solution of eq 17 in conjunction with eq 7 and 15. 
Exact solutions will not be possible in general, but approximate solu- 
tions can always be made. The solutions for <£ will be of the form 
(j)=x±iy and the corresponding values of X will be of the form 
\ k =lj, k ±2iriv k . It should be remarked that the damping constant 
n k is now different for different components k. In general, it will be 
found that lumping the damping resistance in the terminal section 
results in a damping factor for the fundamental and lower harmonics 
greater than would have been obtained had the same resistance been 
uniformly distributed throughout the network. On the other hand, 
the damping factors for the higher harmonics will be found to be much 
less for the lumped damping resistance than for the distributed re- 
sistance. A somewhat more detailed discussion of the solutions of 
such equations will be found following eq 46 where the solution of a 
similar equation is discussed in connection with surge-current circuits. 
As is shown there, the circuit constants in any particular case may 
suggest valid simplifications and approximations which greatly 
facilitate the computations. 

Exact solutions for the amplitude factors P k , see eq 20, cannot be 
obtained in general. With a knowledge of the proper form of solu- 
tion, however, and the correct values of the X's, it will be possible to 
compute, by the ordinary methods of wave analysis, the amplitude 
factors P k necessary to fit a given surge-voltage oscillogram. Even 
this meager information may prove valuable, however. It will en- 
able one, from data obtained on the well-defined crest and tail of a 
voltage wave, to reproduce the sharply rising front of the voltage wave 
which is actually of much importance but is often poorly defined on 
the oscillograms. 



* P. L. Bellaschi, Trans. Am. Inst. Elec. Engrs. 51, 936-945 (1932). 
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III. SURGE-CURRENT GENERATOR 
1. STATEMENT OF THE GENERAL PROBLEM 

The problem to be solved in connection with some types of the 
surge-current generator can be stated conveniently in terms of the 
simplified circuit shown in figure 3. The n identical capacitors com- 



n^nm-w/H-nm — Wr- 



Qz 



Qj 






Qn,. 



Qh 



6 
o 



Figure 3. — Schematic diagram of a surge-current generator. 
The charging circuits have been omitted to avoid confusing the figure. 

prising this generator are connected in parallel by identical leads 
having resistances R and inductances L. The nth, or terminal, 
mesh, which in practice contains the apparatus under test as w^ell 
as other measuring equipment, is closed through a resistance R n 
and an inductance L n which is, in general, different from the resist- 
ance and inductance (JS, L) of the other n-\ meshes. As actually 
used, a surge-current generator will consist of two or more units, 
such as shown in figure 3, connected in parallel to a common sphere 
gap and load, as are shown in figure 4. 6 If the parallel circuits of 
figure 4 are identical, and if we can neglect the effect of mutual 
inductance between leads, it is obvious that the currents in the 
parallel circuits will be identical at all times and that the current 
through the load resistance will be simply three times that flowing 
in any one of the parallel circuits. The solution for the case of a 
symmetrical circuit, such as figure 4, therefore reduces essentially 
to the solution of the circuit of figure 3. As before, we neglect the 
internal structure, if any, of the capacitors and represent their 
internal losses by a suitable series resistance, and neglect the change 
of resistance of the leads with frequency. If identical charges Q 
are placed on each of the capacitors and the nth mesh is then closed 
through its inductance and resistance, it is required to find the 
charges on the vaiious capacitors and the currents in the various 
meshes as explicit functions of time. We are particularly inter- 
ested here in the current in the 7ith mesh. The method of solution 
will closely parallel that used in connection with the surge-voltage 
generator in section II. 



« P. L. Bellaschi, Elec. Eng. 
38, 126(1935). 



3, 86 (1934). K. B. McEachron and J. L. Thomason, Gen. Elec. Rev. 
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2. FORMAL SOLUTION OF THE GENERAL PROBLEM 

Kef erring now to figure 3 and writing down the equations for the 
summation of the voltages around each mesh we have, if Q s is the 
charge on the sth capacitor and i s the corresponding current: 



Q l /C-Ldi 1 /dt-Ri l - Q 2 /C= (T 
Q n /C—L n di n /dt—R n i n =0; 



(26) 



R L R L 




Figure 4. — Schematic diagram of a surge-current generator using three identical 

circuits in parallel. 

Now, noting the currents at point 6, we can write 

ii=— dQi/dt 
i 1 —i 2 —dQ 2 /dt=0 



or 



i 2 =-(dQ 1 /dt+dQ 2 /dt) 



and, in general, 



i s =—2dQ r /dt. 



(27) 
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Let us now make the usual exponential substitution 

Qs=B 8 e-" (28) 

where B s and X are undetermined constants (see footnote 2, p. 588). 

Substituting eq 27 and 28 in the typical equation of eq 26 and 
making a few obvious rearrangements, we have 

(\ 2 LC- \R0)2B r + (\ 2 LC-\RC+ l)B s -B s+1 - 



This suggests the following definitions: 

u=\ 2 LC-\RC+l 1 

Y=\ 2 L n C-\R n C+lj 

Eq 26 can then be put in the form: 

uB 1 -B 2 =0^ 

(u-l)2B T +uB s -B s+1 =0 

r-1 

(F-lJ25 r +m=0 



(29) 



(30) 



The condition for a solution giving values forlthe B's other than 
zero is, as usual, that the determinant of the coefficients shall vanish. 
This nth-order determinant, which leads to the secular equation, we 
shall call F n . The determinantal equation is: 



F = 



u 

u-l 

u-l 


-l 

u 
u-l 


.. 
-1 . . 

u . . 


.. 
. . 
.. 


u-l 

r-i 


u-l 
Y-l 


u-l . . 

r-i .. 


..u-l 
.. Y-l Y 



=0 



nth order. 



(31) 



Let us also define the following determinant which*is obtained on 
expanding F n : 



K n = 



u-l 
u-l 

u-l 


-I 

u 
u-l 


. . 

-1 . . 

u . . 


. . 
.. 
.. 


u-l 
Y-l 


u-l 
Y-l 


u-l . . 
Y-l .. 


..u-l 
. . Y-l Y 



(32) 



nth order. 



Expanding eq 31 on the first row, we obtain 

F n =uF n ^+K n ^ 
Expanding eq 32 on the first row, we obtain 
K n =(u-l)F n ^+K n ^ 
whence K n -i = -F n . 2 + (uF n _ 2 +K n _ 2 ) 

=F n _i—F n _ 2 . 
Therefore, F n =(u+l)F n - % -F n ^. 



(33) 
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Eq'33 is of exactly the same form as eq 14 and the method of solution 
of the two equations is identical. If we make the definition 

0+l)/2 = cosh£ (34) 

we find that eq 33 can be put in the form 

w YsinhnZ— sinh (n— 1)£ faKi 

Fn== slShl (35) 

The problem of solving eq 26 is now reduced to the problem of finding 
those values of £ which will reduce the right-hand member of eq 35 to 
zero. There will be n such values of £ corresponding to the n degrees 
of freedom of the system. 

Values of X are obtained by substituting these values of £ in eq 29 
and 34 from which we obtain 

\ ktl) \ k , 2 =Bl2L±^R 2 /4U+ (4/ZO) sinh 2 ^J (36) 

there being n values of k corresponding to the n possible values of £. 

It can now be verified by direct substitution that eq 30 are satisfied 
for all possible values of £ if the B's are of the form 

B s =U k cosh («-l/2)&. (37) 

The final solution of our problem is therefore given by 

Qs^lU^e-^+U^e-^] cosh (*-l/2)fc (38) 

I 

s 

and i s =— 2 dQJdt 

r = l 

where the U k s are constants to be determined by the given initial 
values of charge and current. 

3. ILLUSTRATIVE CASE, EXACT SOLUTION 

The application of the preceding formulas is best illustrated by 
means of a specific example. For simplicity, a case has been chosen 
in which an exact solution is possible. Other special cases leading 
to exact solutions can be discovered with little difficulty. 

Referring to figure 3, let us assume that 

R n =R and L n =L. 

This is assuming that the terminal section is identical with the typical 
section. Substituting these values in eq 29 and 35, we obtain 

Y=u=2 cosh £—1 
and the condition for a solution, F n =0, reduces to 
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Hence ^-SS+T* 

or, if $=# ^=^+1' (40) 

&=0, 1,2, ...... . Ti-1 

Referring to eq 36 we have 

X,,i,^, 2 -i2/2L±[i? 2 /4i 2 -(4/iC) sin 2 |*T. 

As before, we distinguish three cases depending upon whether R is 
greater than, equal to, or less than the value necessary to critically 
damp a given frequency. For values of R less than the critical 
damping value, we will have 

where n=R/2L (41) 



2 7 rv k =S-R 2 /4L 2 + (4/LO) sin 2 ^? 



If now the value of R is such as to critically damp the rth frequency 
component, we can rewrite eq 38 in a form similar to eq 24. 

Qs= r i2[U k>l e-^+U k ,2e-^] cos (s-1/2)^ 

fc-0 

+ [U T J+U r , 2 ]e-^ cos (s-l/2)* r (42) 

+S C7* cos (s-1/2)^-^ cos (2ttM-^). 

The constants U and /? are determined, as before, by the initial 
conditions, which are, in this case at £=0 

Q s =Qo 

HR is less than the critical damping resistance for the fundamental, 
these conditions lead to the following values for the constants, 

ftt= tan -V/2tt^ (43) 

U k =[(-l) k 2Q /(2n+l) cos ft] cot § 
whence 

&=2(-l)*[2£ /(2n+l) cos ft] cot § cos (s- 1/2) ^<r^ cos (2™**- ft) 



and finally 



3l*/ ™a P.. nr.f2 l±»-uX ci'r 



in=-2Q /(2n+l) (i(7)«S sin ^/ cos ft- cot 2 *&-* sin 2wv k t. (44) 
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Let us assume, as a practical example, a generator having the 
following constants, 

n=10,£=lX10-«A 

C=2.4X10" 7 /,i?=0.01 ohm. 

These constants approximate those that might be used in an actual 
generator. Using these data we find 

M =5X10 3 

V;fc =l^[1.67sin 2 (2A;+l)7r/42-25X10- 6 ] 
Zir 



or 



z^650 sin (2k + l)ir/42 kc/s approx. 
&=0, 1, 2, 9 
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Figure 5. — Current surge from a 10-section surge-current generator. 

Curve I, distributed resistance of 0.01 ohm per section, fundamental very much underdamped. 
Curve II, distributed rasistance of 0.61 ohm per section, fundamental critically damped. 
Curve III, the critically-damped fundamental corresponding to curve 27. 

This gives a fundamental, v , of 48.5 kc/s, which is to be compared 
with 32.4 kc/s, which would be computed by considering all the con- 
stants of the generator lumped and computing the resulting frequency. 
In figure 5, curve /, i n , the current in the terminal mesh, is plotted 
against time, using the circuit constants just given. In figure 5, 
curve II, i n has been plotted for a generator having the same 
constants, except that 5=0.61 ohm. This value of resistance is 
just sufficient to critically damp the fundamental. In this case /* 
(critical) =30.4 X'10 4 . Curve /// represents the critically-damped 
fundamental for this case. 
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4. ILLUSTRATIVE CASE, APPROXIMATE SOLUTION 

In actual practice it is customary to make the internal resistance of 
the generator (R) as small as possible and to achieve the necessary 
damping by making the terminal resistance (R n ) sufficiently large. 
These conditions allow certain simplifying assumptions to be made in 
the general solution of eq 35. It is well known, and can be verified in 
the numerical example of the last section, that if the damping factor 
is small, the damping affects the frequency as a second-order correc- 
tion only. Neglecting the internal resistance of the generator for the 
moment, we have, from eq 29 and 34, 

X=±2/(Z<7) 1/2 .sinh£/2 (45) 

and eq 35 then becomes 

cosh (n-l/2)f+[2£ n /Z.sinh £/2Tfl„. (<?/£) 1/2 ] sinh n£ n fAAX 
c^shl/2 =0 (46) 

This equation can be solved approximately by any one of the standard 
methods. Its roots will be of the form 

and the X's will be of the form 

X..i,X., 2 -±2/(i(7) 1/2 -sinh^i-^±27r^, 

and we have, in addition, the small damping arising from the inter- 
nal resistance /z — R/2L. 

If we make the additional assumption that L n =L, that is, that we 
add pure resistance only to the terminal section, eq 46 simplifies to 

cosh (n+l/2)ZTR n (C/Ly< 2 sinh nj 

cosh £/2 ^ ' 

Approximate solutions can be written down for two limiting cases 

2k4-l 
R n very small ^=K~TA ir ^ = ^ 1 1? 2, . . . n—l 

k 
R n very large f*=-xi k=l, 2, .... n—l 
n> 

Taking one of these approximations as a value of iy above, we can 
enter eq 47 and find approximate values for x k . Using these values of 
x k , we can now compute a second approximation to the true values of 
y k . This process may be repeated if necessary. The value of y cor- 
responding to k— may sometimes have to be obtained by a separate 
evaluation. If we wish to find the value of R n necessary to critically 
damp the fundamental, we solve equation 47 graphically under the 
assumption that y is a small quantity the squares of which can be 
neglected and find that value of i?i which just causes y to vanish. 

89799—36 8 
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Sample calculations have been carried out on the 10-section 
generator of the last section under the conditions L n =L, I? =0.01 
ohms, B n =2.3 ohms (critical damping for the fundamental). In 
figure 6 are plotted, as solid lines, the damping factors, p, k} for this 
case against the corresponding frequencies. As dotted lines are also 
plotted, the corresponding damping factors against frequency for the 
case of the same generator, the fundamental of which is critically 
damped by distributed resistance, R=R n =0.Ql ohm. 

It should be remarked that had we considered the constants of this 
generator to be lumped, we would have computed a terminal resistance 
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-Damping factors against frequencies for a 10-section surge-current 
generator. 



Solid lines, damping factors against frequencies for the case of a lumped terminal resistance just sufficient 
to critically damp the fundamental. i?=0.01 ohm, ^=2.3 ohms. 

Dotted lines, damping factors against frequencies for the case of a distributed resistance just sufficient 
to.critically damp the fundamental. R=R n =0M ohm. 

of 4.1 ohms necessary to critically damp the system. The damping 
resistance obtained from the present closer calculations, that is 2.3 
ohms, is in reasonable agreement with the experimental values 
observed on just such a generator, that is 1.5 to 2.0 ohms. 7 

IV. CONCLUSION 

It is thus seen that the mathematical theory of surge-generators 
presents no essential difficulties, and that it is possible to set up com- 
plete formal solutions for both the surge- voltage and the surge-current 
generator which take account both of the internal structure of the 
generator and of the load to which it is connected. From these solu- 
tions it is possible to obtain exact numerical values for the damping 
factors and frequencies in certain special cases, and to obtain approxi- 
mate values in all cases. Exact expressions for the amplitudes and 
phase angles can be obtained only in certain special cases. 

» P. L. Bellaschi, Elec. Eng. 53, 86 (1934). 
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It has been the purpose of this discussion to outline the general 
method of attack on the mathematical problem presented by surge- 
generator circuits and the difficulties which will be met rather than to 
discuss in detail all the pecial cases which might be derived from the 
equations given above. If a particular generator were under con- 
sideration, many valid assumptions and approximations leading to 
useful approximate solutions might suggest themselves, which it 
would be fruitless to consider in a general discussion. 

Washington, July 23, 1938. 



